Surface spectral function in the superconducting state of a topological insulator 



O 

(N 

oo 



o 



> 
\o 

o 



X 



Lei Hao and T. K. Lee 

Institute of Physics, Academia Sinica, NanKang, Taipei 11529, Taiwan 
(Dated: January 13, 2013) 

We discuss the surface spectral function of superconductors realized from a topological insulator, 
such as the copper-intercalated Bi2Se3. These functions are calculated by projecting bulk states to 
the surface for two different models proposed previously for the topological insulator. Dependence of 
the surface spectra on the symmetry of the bulk pairing order parameter is discussed with particular 
emphasis on the odd-parity pairing. Exotic spectra like an Andreev bound state connected to the 
topological surface states are presented. 
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I. INTRODUCTION 

Recently, a new state of matter, the topological 
insulator, attracts much theoretical and experimental 
attention. [H, i i, i, i i, 0, i, i M, El, El, Hi, HI, P 
It has an electronic structure dominated by the spin- 
orbit coupling, which is a band insulator with a well- 
defined gap in the bulk but can host an odd number 
of Dirac cones protected by time reversal symmetry on 
the surface. [l|, H, 0| The intrinsic spin-orbit coupling 
makes it promising for spintronics applications. [lOj The 
proximity induced superconducting state on the topo- 
logical insulator surface by a deposited superconductor 
was proposed to create Majorana fermions^lGj], which 
may provide a new wa y to realize the topological quan- 
tum computation. [13, llSl j Lately, the realization of a 
superconducting state in a typical topological insula- 
tor Bi2Se3 by intercalating Cu between adjacent quin- 
tuple units_(Cu^Bi2Se3) makes the system even more 
attractive. [1^ [2^ The large diamagnetic response shows 
that the pairing is mainly of bulk character. In another 
topological insulator Bi2Te3, the application of a high 
pressure also turns the material into a superconducting 
state, [m, HI] 

Since the superconductivity is bulk and intrinsic to the 
material, if zero energy surface Majorana fermion mode 
exists, it would be easier to manipulate as compared to 
that induced by proximity effect, as it would not expe- 
rience the interface roughness or mismatch common to 
a junction type device. Possible nontrivial odd-parity 
pairing in Cu2:Bi2Se3 is proposed and analyzed by Fu 
and Berg. It was argued that only if a bulk gap opens 
and the bulk pairing is odd in parity, would zero en- 
ergy Andreev bound states appear in the surface spec- 
trum. However, their analysis was concentrated on the 
case when the chemical potential is much larger than the 
gap and the topological surface states are already merged 
into the continuum conduction band. A similar analysis 
was put forward by Sato. [24] 

Despite the above works, a detailed theoretical analysis 
of surface spectrum in the superconducting phase arising 
from a topological insulator is still lacking. In particular, 
not much is reported for the situation when the chem- 
ical potential is only slightly larger than the insulating 



gap and both topological surface states (or, the surface 
conduction band[20l]) and the continuum bulk conduc- 
tion band are present but separated. Since the surface 
spectrum is central to the topological properties of a ma- 
terial both in the normal and in the superconducting 
state, which is also directly accessible by experimental 
techniques as ARPES[20] and STM 25], it is highly de- 
sirable to make a detailed study of them. This will help 
to understand better superconductivity in systems with 
nontrivial topological band structure. 



We focus on two questions in this paper. One is the 
effect of superconducting pairing on the topological sur- 
face states [20] present in the normal state. The other is 
the existence of surface Andreev bound states. We have 
noticed that two different models 12, 26] are often used 
in-discriminatively in literature. They have the same nor- 
mal state energy spectra but may be different in the su- 
perconducting state. We thus present our results for both 
of them. What happens to the topological surface states 
when pairing is introduced in the bulk depends on the 
orbital character of the topological surface state and on 
the bulk pairing symmetry. Only when the continuum 
part of the band opens a full gap and the topological 
surface states, while separated from the bulk conduction 
band, do not open a gap for an odd-parity pairing, would 
a gapless Andreev bound state appear. We find that the 
orbital characters of the topological surface states are dif- 
ferent for the two models. For a certain bulk pairing sym- 
metry, it is possible that the topological surface states of 
one model opens a gap while that of the other model is 
still intact. The existence of Andreev bound states, the 
most important indication of nontrivial topological order 
in the superconducting phase, is thus expected to be also 
related to the orbital character of the topological surface 
states. We show that the interplay between the contin- 
uum bulk conduction band and the topological surface 
states produces a ring or a segment of zero energy states 
in addition to the Andreev bound states depending on 
the symmetry of bulk pairing order parameters. 
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II. MODELS AND THE NORMAL STATE 
SURFACE MODES 

In the fohowing, when we talk about topological in- 
sulators, we would be mainly referring to Bi2Se3, which 
shows a well defined Dirac cone structure for the topo- 
logical surface states. [ill, [13, US] The model we consider 
below could be easily generalized to study other topolog- 
ical insulators like TlBiSe2 [H, [ll] and Bi2Te2Se[3a [3l[. 

The band structure of Cua;Bi2Se3 is similar to that 
of Bi2Se3, the most essential part of which consists of 
two Pz orbitals on the top and bottom Se layers hy- 
bridized with neighboring Bi pz orbitals, in each quin- 
tuple Bi2Se3 unit. ,12. , 23,] In the presence of spin-orbit 
coupling, the normal state has four degrees of freedom. 
Label the two orbitals concentrating mainly on the top 
and bottom (seeing along the —z direction) Se layer of 
the various Bi2Se3 quintuple units as the first and second 
orbital, the basis is taken as i/'k=[cikt, C2kt, cik^, C2k^]"^- 
The models could be written compactly in the following 
matrix formlll, [H lH 



(1) 



i7(k) = eo(k)/4x4 + "^»(k)r. 



I4x4 is the fourth order unit matrix, giving rise to a topo- 
logically trivial shift of the energy bands and would be 
neglected in most of our following analysis. In terms of 
the two by two Pauli matrices Si (i=0, • • • , 3) in the 
spin subspace and (Ji (i=0, • • • , 3) in the orbital sub- 
space, the first three Dirac matrices are defined as[T2il26j 
ro=So(8)cri, Vi—si®a^, T2—S2®cr^. As regards Fq, there 
are presently two different choices: (I)so cr2 2^ [1^ HI] 
and (II) S3 (g) crafU H, m [M [li, which define the two 
models that would be considered in parallel in the fol- 
lowing discussion. Since Bi2Se3 is inversion symmetric, 
parity could be used to label states. Some papers adopt 
the bonding and antibonding states of the two orbitals 
defined above as the orbital basis. [H, HI, [13, [13, [IE [IB] 
The corresponding models could be obtained in terms of 
a simple unitary transformation performed in the orbital 
subspace. Note that, the two different models give the 

same bulk band dispersion e± (k) =eo (k) ± \J Yl^i=o C*) ' 
with each of two eigen- energies two fold Kramers degen- 
erate due to the time reversal symmetry and the inversion 
symmetry of the model. However, we will show in the 
following that they are essentially two physically distinct 
models. 

For the coefficients eo(k) and TOi(k) [i — 0, 1, 2, 3), 
there are different possible parameterizations which co- 
incide with each other close to the F point. [T2l. [T^ [26l W\ 
Without loss of generality, we take the parameterizations 
of Wang et al. [26| . Since the diagonal term eo (k) propor- 
tional to the unit matrix does not affect the topological 
characters, it is ignored here. The system is defined on 
a hypothetical bilayer hexagonal lattice stacked along 
the z-axis, respecting the in plane hexagonal symmetry 



of the original Bi2Se3 lattice. With the three inde- 
pendent in-plane nearest neighbor unit vectors defined 



as 6i=(^, 



(-^, 5), and 63^(0, -1), we have 

mo(k)=m + 2tz{l — cos kz) + 2t{3 — 2 cos ^k^ cos ^ky — 

cosky), mi(k)=2\/3t sin ^fcj, cos j/cj^, 

m2(k)— 2t{cos^ kx sin ^ky + sinfcj,), and 
m3(k)=2tz sinkz- The in-plane and out-of-plane 
lattice parameters [ssj are taken as length units in 
the above expression, that is a=c=l. When mtz<0 
and mt<0, the parametrization defined above [26i and 
the parametrization in the small k effective model 
proposed by Zhang et a/. [12] describe qualitatively the 
same physics. With this parametrization, it is easy 
to see that the model has the inversion symmetry 
PH{'k)P~^=H{—'k), where the inversion operator is 
defined as P=so ® cri.[l] 

Now, we clarify the differences between models (I) and 
(II) by their surface states, which is one of the most 
important signatures of nontrivial topological order in 
the system. Close to the F point in the BZ, we take 
TOi=o, - .3(k)={TO + ^t{kl + kf) + tzkl^itkx, 5tky, 2tzkz}, 
in which t > 0, tz > and m < 0. Consider a sample 
occupying the lower half space z < 0. The possible sur- 
face states localized close to z = is searched by solving 
a set of four coupled second order differential equations 



Hikx 



Q,kz 



-idz)-^{z) = S*(z), (2) 



together with the open boundary condition 

^{z)\z^o=^{z)\z=-oo=Q.m. [33 *(^) is the four- 
component eigenvector and E is the energy of the surface 
mode, respectively. We look for the zero energy states 
and hence set i?=O.0] 

For the model (I) with F3=so®o'2, the up and 
down spin degrees of freedom are decoupled from each 
other. The wave function could thus be written as 
^'(z)=Jui(z), M2(z), U'i{z),Ui{z)]^^[x^{z),Xi{z)Y . The 
two spin components of the zero energy mode satisfy the 
same equation as (s is f or J, for the two spin degrees of 
freedom) 



l(m - tzdl)ai - 2ttzdz<T2]xs(z) = 0. 



(3) 



The two degenerate zero energy surface modes for z<0 
are obtained as 



*a(;2) = C77„(e^/«+ -e^/^-). 



(4) 



where a=l or 2, C is a normalization constant and 
iZ^=l ± y/l + m/tz- 1/Re[^±^] ('Re' means taking the 
real part of a number) are the two penetration depths of 
the surface modes into the bulk. The two unit vectors 
are {r]i)p=S^i and {112)13=6^3, where Jq.^ is one for q;=/3 
and zero otherwise. Take {"i/i, ^'2} as the two basis, the 
effective model for the surface states are obtained by con- 
sidering the kx and kj^ dependent terms in the original 
model as perturbations, which are 



AH: 



3D 



-Akl 



kl)To 



3tikxri + kyT2). (5) 
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Suppose the two basis are normalized, the effective model 
for the surface states isfl^l 

HeffCk) = it{k.r,s^ + kySy), (6) 

where Sx and Sy are the first and second Pauli matrices. 
Since the two basis both have definite spin characters, Sx 
and Sy in the above equation could also be considered as 
acting in the spin subspace. The most salient feature 
of this model is that the corresponding surface states 
has contributions only from the first orbital. When we 
consider a sample occupying z>0, the surface states at 
z=0 would arise only from the second orbital. 

We now study model (II) for r3=S3(g)CT3. We still con- 
sider a sample situated at z<0 with the open boundary 
conditions. Following exactly the same steps as for the 
first model, we obtain the two degenerate zero energy 
surface states as 

<fa = Cr,„(e^/«+-e^/«-), (7) 

where a—1 or 2, C is a normalization constant and ^± 
are defined identically as above. However, the two unit 
basis vectors are quite different from the first model and 
are ?7i = ^[l, -i, 0, 0]"^ and r/2=^[0, 0, -i, 1]^. The two 
dimensional effective model for the surface states is also 
a bit different at least formally, which is 

Heff{k) = 3tz • (k X s) = 3t{kxSy - kySx). (8) 

The Pauli matrices Sx and Sy act in the two fold degen- 
erate basis of the zero energy surface states which both 
have definite spin characters. For a general two dimen- 
sional wave vector, the surface state would be a linear 
combination of all the four spin-orbital basis. 

Thus there are qualitative differences between two 
models which were used in-discriminatively in the litera- 
ture for BiaSea . [13, H, IM IH, H, HI HI, While only 
one orbital contributes to the surface states for model 
(I), both two orbitals contribute in equal weight to the 
surface states for model (II). On the other hand, the effec- 
tive model of the surface states has the same spin-orbital 
coupled form as the linear kx and ky terms in the origi- 
nal three dimensional model for model (I). However, the 
effective model is changed from k • s to ^ • (k x s) for 
model (II). We have verified that, if we change the in- 
plane spin-orbit coupling of the two models (I) and (II) 
from k - s to z • (k X s), the resulting effective model of the 
surface states would have the form z • (k x s) for model 
(I) but will be switched to k • s for model (II) . 

Before ending this section, we would like to point out 
that both the fc^-linear term in (k) and the fc^-square 
term in r7io(k) are essential to obtain the zero energy 
surface modes. If we omit the k1 term in mo(k), it is 
easy to verify that the gapless surface states no longer 
exist. This is related to the fact that band inversion 
is essential to the appearance of nontrivial topological 
surface states [3,[3^, which can only occur in the presence 
of k1 term for kx=ky=Q. 



III. SUPERCONDUCTING STATE SPECTRAL 
FUNCTION 

A. surface Green's functions in the 
superconducting state 

The realization of the superconducting state in 
Cu:EBi2Se3[ll, HO] has brought about excitement that 
non-trivial topological superconducting state might be 
realized in this system, in which topologically protected 
gapless surface states traverse the bulk superconducting 
gap.jH, [23, The recent realization of superconduct- 
ing phase in Bi2Te3 under high pressure [21I [23 | makes the 
Bi2X3 (X is Se or Te) material a very promising candidate 
system to realize topologically nontrivial superconduct- 
ing phases.HHIiO] 

The normal state of the topological insulator is marked 
by the presence of topological surface states inside the 
bulk gap. These gapless topological surface states are 
well separated from the bulk conduction band at low en- 
ergies and become indistinguishable for energies much 
higher than the conduction band minimum. Depending 
on the doped charge density the superconducting state 
could occur with the chemical potential either deep in the 
bulk conduction band or in the intermediate region where 
the topological surface states are well separated from the 
bulk conduction band'2^. In the latter case, the cou- 
pling between the continuum bulk states and the isolated 
topological surface states may cause some new interesting 
phenomena. Thus this intermediate region is where we 
will concentrate on below. Furthermore, the actual pair- 
ing symmetries of the superconductin g C u3,Bi2Se3 and 
Bi2Te3 are presently unknown. [H [Hill ^ Thus we 
will examine cases with different pairing symmetries in 
the hope to provide clues to identify the pairing symme- 
try and the role contributed by the topological surface 
states. 

In the following, we will study the surface spectral 
function to see possible nontrivial topological properties 
arising from the normal phase topological order which 
is subject to a certain bulk pairing. The surface spec- 
tral function, which could be obtained from the surface 
Green's function, has been studied by ARPES[23 and 
STM^l to give important information on the topological 
properties of the system. In the superconducting state, 
we expect to see some surface Andreev bound states if a 
certain superconducting order is realized in the material. 

In the presence of a surface perpendicular to the z-axis, 
kx and ky are good quantum numbers, and k^ is replaced 
by —idz as we shall search for surface states. We then dis- 
cretize the z coordinate and turn the whole sample (z<0) 
with a surface at z=0 to a coupled quintuple-layer sys- 
tem. Label each separate quintuple unit with an integer 
index n, and make the substitutions dz'4'n{z) — \[4'n+i ~ 
i/'n-i] and 9f'0„(z)=-0n+i + i'n-i - (c is set as 
length unit along z axis), the Hamiltonian consists now 
the intra-lajer terms and the interlayer hopping terms, 
H=H\\ + H±. The intra-layer part of the model is 
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-^ll=I]«k^nk'*a;2/(k)V'„k, in which 

h,y{k) - TOj,(k)ro + 77ii(k)ri + m2(k)r2. (9) 

r7ii(k) and r7i2(k) are the same as those in the bulk 
model. ?TiQ(k) is obtained from r7io(k) by first expanding 
it up to the square term of kz and then replacing the 
term proportional to k'^ from Ck"^ to 2C. 

The inter-layer hopping term is 
-^-L=X]nk V'lk'*zV'n+i,k+H.c. In tcrms of the pa- 
rameterizations of Wang et al. [1^ , we have 



Since we have now two orbital and two spin degrees 
of freedom, there are many possible pairing channels for 
different possible pairing mechanisms. Realistic theoret- 
ical determination of the pairing symmetry requires the 
knowledge of pairing mechanism and reasonable param- 
eter values, which are both lacking presently. [2j| Here 
we will consider singlet and triplet pairing orders as phe- 
nomenological input parameters. Their qualitative dif- 
ferences in spectral functions could help to identify the 
pairing symmetry from experiments. 



hz = ~tz{TQ + iVz 



(10) 



Since the matrix now appears only in the hz part of 
the coupled layers system, the difference between the two 
models enters only through the interlayer hopping term. 

Now introduce superconducting pairing and define the 
Nambu basis as ^Ik— I^nk' V'n-kl- The intra-layer part 
of the Bogoliubov de Gennes (BdG) Hamiltonian is then 



sc_ 



^i?sc(k)(/)„k, hi which;41] 



hoik) A(k) 
-A*(-k) -h*{~k) 



(11) 



where hQ{'k)=hxy{k) — fil^xi, with fi the chemical po- 
tential. A(k) is the 4x4 pairing matrix. Ignoring the 
possibility of interlayer pairing, the interlayer hopping 



terms are Hf^- 



At 



Hz(t)n+i.k+ii-C-, in which 



Hz = 





-h* 



(12) 



Once the pairing order is given, the surface spectral 
function is obtained from the retarded surface Green's 
functions, which could be calculated in terms of stan- 
dard transfer matrix method. (43| In the simplest form 
of the method, the 8x8 retarded surface Green's func- 
tion G(k, Lu) is obtained by self-consistent calculation of 
G(k, w) and a transfer matrix T(k, w) as[26l [i^ 
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hIt. 



T = GHz 



(13a) 



(13b) 



where g=[z/8x8 — Hsc{^)]~^ {z=oj + i-q) is the retarded 
Green's function for an isolated layer, r/ is the positive 
infinitesimal 0^, which is replaced by a small positive 
number in realistic calculations. Self-consistent calcula- 
tion of the Green's function starts with G=g. The surface 
Green's function could also be obtained in terms of other 
iteration schemes, such as the algorithm in Rcf.[42]. We 
have found no difference between the results obtained in 
terms of different iteration schemes. 

After the retarded Green's functions are at hand, the 
spectral function is obtained as 
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A(k, w) = - 2^ ImGii (k, i:j)/7r. 



(14) 



B. gap opening in the topological surface states 

Before presenting the full spectral function, we first 
would like to examine what happens to the topological 
surface states inherited from the normal state _20i upon 
the formation of a certain bulk pairing. The most salient 
feature of the topological insulator is the presence of gap- 
less surface states (3D) or edge states (2D).P, |^, 0, HI] 
In the case of Cua;Bi2Se3, it is found that these surface 
states in the non-superconducting Bi2Se3 persist to the 
superconducting copper intercalated samples and are well 
separated from the bulk conduction band and hence well- 
defined. [19|, i20|] It is thus an interesting question what 
would happen to them if a certain pairing forms in the 
bulk. In this subsection, we give a simple criterion to 
judge whether a gap would be induced in the topological 
surface states for an arbitrary bulk pairing. 

Suppose the chemical potential lies slightly above the 
bottom of the bulk conduction band where the topologi- 
cal surface states are well separated and well defined. [2^ 
If a pairing is realized in the topological surface states, it 
should occur between the two time reversal related states 
for k and -k.fl?'] 

We first consider model (I). For our purpose, we would 
concentrate on the positive energy branch of the topo- 
logical surface states. When the pairing occurs in the 
valence bandj21„ i23], the analysis and conclusion would 
be similar. Since pairing occurs in the (k^:, kj,) space, 
we would ignore the z-dependence of the surface modes 
when analyzing pairing properties. From the basis and 
the effective model obtained in Sec. II, the two eigenvec- 
tors for a certain 2D wave vector are 



7y„(k) = -L[l,0,a^,0]^ 



(15) 



where a is '-I-' ('— ') for the upper (lower) branch of the 



surface states, k±=k^ zLiky, k—^k'^ 
tion operators of these states are 



ky. The annihila- 



V2 



[cikt ■ 



ak^ 



(16) 



If pairing is induced in the upper surface conduction band 
at k, the only possible pairing would be proportional to 
d|^^(i^j,^.[17] Denote the time-reversal operator as T.fT?!. 
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reversal symmetry of the pairing, the actual pairing 
should be of the form 



m Since r4ktT-i=c|_ki and Tc\^^^T-^ 
have Tdt,d^_^,T^^ =i^dt , d\, , . To ensure the time- 



Ao'fdl^dl 
2 ^" 



k ikt'-i- 



+(cl 



k^"^l-kt 



_ _ t t 
-kt ^ "^Iki"^!- 

^Ikt'^l-ki)]' 



■ki 



(17) 



where Aq is the real pairing amplitude, which could be an 
even or odd real function of k depending on the pairing 
realized in the bulk. 'SCB' is abbreviation for the surface 
conduction band (the topological surface states). Thus, 
the surface conduction band only supports the anti-phase 
Px ± ipy equal-spin triplet pairing and the spin-singlet 
pairing within orbital 1. No other bulk pairing channels, 
especially those inter-orbital pairings, would open a gap 
in the topological surface states within the framework of 
model (I). 

For model (II), the two eigenvectors of the surface 
states for a certain 2D wave vector are (again, ignoring 
the z-dependence) 

,y„(k) = i[l,-^,a^,^a^]^, (18) 

where a is ('— ') for the upper (lower) branch of the 
topological surface states. Following the same arguments 
as for the first model, when the chemical potential cuts 
the upper branch of these well defined surface states the 
time reversal invariant pairing is of the form 



^ I ^ I'-lkt'^l-kt ''2kt^2-kt'' 



-^(c;,,,c 



t J 



^2-kl^ 



. k^ / t t t t \ 

+*~^i'^lkt"^2-kt + '^2kt'^l-kt) 

■ ('ft t t ^ 

+*~^(.'^lki'^2-k^ + '^2k^'^l-k^) 

+ ("^lk^'^l-kt ~ '^Ikt'^l-ki 
+'^2k^''2-kt ^ 4kt"^2-ki) 



lkt^2-k^ ' ^ii^i'-2-\i'[ 

,t J J J 



^2kt^l-k^ '^2k^^l-kt 



"l-kt)]- 



(19) 



As in Eq. (17), Aq could be a constant or a real function 
of k compatible with the symmetry of one pairing com- 
ponent contained in the above decomposition. Besides 
the intra-orbital pairing channels active in the model (I) , 
there are two additional inter-orbital pairing channels 
that are effective in producing a gap in the topological 
surface states. The last term in the above equation is 



just the odd-parity inter-orbital triplet pairing proposed 
by Fu and Berg [23] as a possible candidate of a topolog- 
ical superconductor to be realized in a superconductor 
like Cu:EBi2Se3. 

The clear difference between Ag^^ and A^^^ makes 
the distinction between model (I) and model (II) more 
obvious. Since the gap opening of the surface states 
is measurable, it is highly desirable to ascertain which 
model is the correct description of the underlying physics 
of Bi2Se3 and Bi2Te3. 

Previously, a simple effective model calculation indi- 
cates that no gap opens in the topological surface states 
for any triplet pairing induced by proximity effect on the 
surface of a topological insulator. 4l| However, our anal- 
ysis above indicates that if the proximity induced triplet 
pairing is compatible with any of the triplet components 
explicit in A;^^^ (A^^^) for model I (model II), then a 
full pairing gap could still be opened in the topological 
surface states. Note that the real gap opening pattern in 
the topological surface states also depends on Aq. 

Except for the pairing channels explicit in A;^^^ for 
model (I) and A;^^^ for model (II), no other bulk pairing 
could open a gap in the topological surface states. The 
existence of surface Andreev bound states depends on 
whether or not a gap opens in the topological surface 
states. We clarify this matter in the next section. 



C. spectral function for typical pairing symmetries 

Observation of superconductivity in Cu2:Bi2Se3 brings 
about anticipation that nontrivial topological supercon- 
ducting states might be realized in this material. The 
topological superconductor is defined as a state with a 
full pairing gap in the bulk and nontrivial gapless An- 
dreev bound states on the surface. [23j 

Possible pairings realizable in a system depend on the 
symmetry of the system and the specific pairing mech- 
anism. In the case of pairing induced by short range 
electron density-density interactions, Fu and Berg iden- 
tified four possible pairing channels. [23| However, if the 
pairing is induced by more long-range interactions, such 
as the electron-phonon interaction, other pairing chan- 
nels (e.g., in which the pairing potential is k dependent) 
would also be possible. In the following we would ana- 
lyze several typical pairings and compare results of the 
two different models. In each case, there are three typi- 
cal situations as regards to the position of the chemical 
potential ^: (1) fj, lies in the bulk gap; (2) ^ lies above 
but close to the bottom of the bulk conduction band, 
where the topological surface states are well separated 
from the continuum bulk conduction band; (3) /i lies far 
above the bulk conduction band bottom, where the sur- 
face states have merged into the continuum conduction 
band. While the latter two cases are relevant to the su- 
perconducting state of Cu2,Bi2Se3[ll, [131, the first case 
could be regarded as mimicking the proximity effect from 
an external superconductor. [TtI [l8l.l4l| In this paper we 
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would focus on the latter two situations. When the chem- 
ical potential lies in the valence band[2TI. [2^. the results 
should be qualitatively similar for the same type of bulk 
pairing. 

Follow Wang et al. [2^ , the model parameters are taken 
as t=t;,=0.5, m=-0.7 in most cases. 0.7 is half of the bulk 
band gap. The width of the bulk conduction band at 
kx=ky=0 is 2{2tz — |to|)=0.6. The small positive number 
ij in the Green's functions is taken as 10^''. 

even-parity intra- orbital singlet pairing First, 
we study the simplest possible pairing denoted by 
A(k)=jAoS2 (8i (To. Spectral functions for the two dif- 
ferent models are the same for this pairing, so only one 
is presented in Fig. 1. Here and in the following, the de- 
gree of darkness indicates the intensity of the spectrum. 
The continuum portions of spectrum are contributions 
from the bulk states, which have small finite amplitudes 
on the surface. Henceforth, they would be called bulk 
conduction band for simplicity. The contributions from 
the topological surface states are somewhat speckled be- 
cause we have taken a finite grid in the (k, oj) plane to 
calculate the spectral function. When the grid points are 
taken to be very dense, contributions from the topological 
surface states will also become smooth. To see the qual- 
itative behavior more clearly, a reasonably large pairing 
amplitude Ao=0.1 is considered. [20j The result is nearly 
identical in the FK direction (along axis) and the FM 
direction (along ky axis) of the 2D reduced Brillouin zone 
(BZ). The topological surface states of both two models 
open a gap, which are consistent with the analysis of the 
previous subsection. Since no Andreev bound state ex- 
ists, this pairing is topologically trivial. The other intra- 
orbital singlet pairings with a k-dependent Aq, which is 
an even function of k could also be considered, such as 



the 



-wave pairing. In these cases, there would be 



line nodes along the nodal directions of the pairing gap. 

- even-parity inter-orbital singlet pairing Since there 
arc now two orbits, another singlet pairing exists in the 
inter-orbital channel. The pairing matrix for the s-wave 
case is A(k)=iAoS2 (8i ui. The corresponding spectral 
functions presented in Fig. 2 for this pairing are still 
identical for the two models. They differ from the spectra 
of the former intra-orbital pairing channel in at least two 
aspects. First, no gap opens in the topological surface 
states, which is in agreement with the criterion proposed 
in the previous subsection. Second, though a full gap 
also opens in the continuum part of the spectrum, it is 
not constant and shows some k dependence. As shown 
in Fig. 2(b), the continuum part of the spectrum even 
nearly closes at some special wave vectors for certain pa- 
rameters. Another interesting feature is the strong redis- 
tribution of spectral weight between the continuum con- 
duction band and the topological surface states. Some 
weight in the bulk conduction band part of the surface 
spectrum above the chemical potential is depleted and 
transferred to the topological surface states below the 
chemical potential. This redistribution arises from the 
particle hole mixing induced by the presence of bulk su- 




FIG. 1: Spectral function for even-parity intra-orbital s-wave 
pairing, for two typical parameter sets for which the topolog- 
ical surface states at the chemical potential (a) are well sepa- 
rated from the bulk conduction band and (b) merges into the 
bulk conduction band. The two models give identical results 
for this pairing. Spectrum along other directions (crossing 
the r point) are qualitatively identical. 



perconducting pairing. As would see below, a similar fea- 
ture is present for each bulk pairing that does not open 
a gap in the topological surface states. 

- odd-parity inter-orbital triplet pairing We now con- 
sider the odd-parity inter-orbital triplet pairing chan- 
nel, proposed by Fu and Berg as a candidate for pos- 
sible nontrivial topological superconducting states in 
Cua;Bi2Se3. 23] The pairing matrix is A(k)=AoSi (8) cr2. 
As was shown in Fig. 3, the spectral functions for the 
two models differ greatly. When the chemical potential is 
close to bottom of the bulk conduction band, the surface 
conduction band is still gapless for model (I) but opens 
a gap for model (H), in agreement with the analysis in 
the above subsection. Another essential difference is the 
existence or not of Andreev bound states. For model (I), 
a band of Andreev bound states appears inside of the 
insulating gap of the continuum which continuously con- 
nects to the topological surface states. While for model 
(H), there is a point node at (0, 0), no Andreev bound 
state exists inside of the gap region. When the chemical 
potential is increased to the position where the surface 
conduction band has almost merged into the continuum 
part of the surface spectrum (corresponding to contri- 
bution from the bulk conduction band), surface spectra 
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0.50 



FIG. 2: Spectral function for even-parity inter-orbital s-wave 
pairing, for three sets of parameters for which at the chemical 
potential (a) the topological surface states are well separated 
from the bulk conduction band, (b) the topological surface 
states are almost merged into the bulk conduction band, and 
(c) the topological surface states are well merged into the bulk 
conduction band. The two models give identical results for 
this pairing. Spectrum along other directions (crossing the F 
point) are qualitatively identical. 



for the two models are as shown in Fig. 4. Since now 
there is no well separated surface conduction band, a full 
gap opens also for the model (I). However, a band of 
Andreev bound states still exists. When we further in- 
crease the chemical potential to /i > 1.3 (for m—-0.7), 
there is no state close to the F point, then there would 
be no Andreev bound state even for model (I). A strong 
redistribution of spectral weight arising from the parti- 
cle hole mixing between the continuum bulk conduction 
band and the topological surface states is observed in 




o.o 



0.1 



0.2 



0.3 



O.4. 



FIG. 3: Spectral function for odd-parity inter-orbital triplet 
pairing in the opposite spin pairing channel, for (a) model 
(I) and (b) model (II). As shown are the cases for which the 
topological surface states are well separated from the bulk 
conduction band at the chemical potential. The parameters 
are as shown on the figures. Spectrum along other directions 
(crossing the F point) are qualitatively identical. 



results for model (I). 

Besides the pairing studied above, there are also other 
odd-parity pairings. Since they are possibly related with 
nontrivial topological superconducting phases, we ana- 
lyze several of them in the following. 

- odd-parity intra- orbital singlet pairing In this case, 
every orbit pairs into a spin-singlet, but the two or- 
bits has a relative tt phase difference and is hence 
odd in parity. [13 The pairing matrix is denoted as 
A(k)=iAoS2 ^ era. The corresponding spectral functions 
for the two models are presented in Figs. 5(a) and 5(b) 
for chemical potential close to the bottom of the bulk 
conduction band and thus the topological surface state 
is well defined. A comparison with the previous pairing, 
shown in Figs. 3 and 4, indicates that the results are 
interchanged between the two models. Gap opening in 
the topological surface states again follows the expecta- 
tion from the previous subsection. The Andreev bound 
states in the bulk band gap for the second model again 
connect continuously to the protected topological surface 
states. The behaviors for higher chemical potentials, as 
shown in Figs. 5(c) and 5(d), are similar to that of the 
previous pairing shown in Fig. 4 with the two models 
interchanged. 
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FIG. 4: Spectral function for odd-parity inter-orbital triplet 
pairing in the opposite spin pairing channel, for (a) model 
(I) and (b) model (II). As shown are the cases for which the 
topological surface states are merged into the bulk conduction 
band at the chemical potential. The parameters are as shown 
on the figures. Spectrum along other directions (crossing the 
r point) are qualitatively identical. 



- odd-parity inter- orbital equal-spin triplet pairing An- 
other interesting possibility is the equal-spin pairing 
channel. Here we consider the two fold degenerate 
inter-orbital pairing channels as proposed by Fu and 
Berg.[23| This kind of pairing could be favored by in- 
terorbital ferromagnetic Heisenberg interactions, [i^ The 
two independent choices for the pairing matrix are 
A(^)(k)=iAoSo <S) (72 and A(^^(k)=AoS3 ® 02- For these 
pairings, it is easy to see that no gap opens in the topo- 
logical surface states for both models. Since results for 
the two models are identical, we only show those for the 
model (I). As shown in Fig. 6 for A^^''(k), a peculiar 
anisotropic Andreev bound state structure is observed. 
An important difference of this pairing from the above 
odd-parity pairing channels is that it is anisotropic with 
respect to kx and ky. Though the bulk dispersion is gap- 
less in the kyk^, {k^k^) plane for A(^'(k) (A'-'^^k)), there 
is still a band of Andreev bound states for the wave vec- 
tors smaller than kp where a gap opens. The peculiar 
feature of the Andreev bound states along ky (k^) for 
A'-^-*(k) (A^^-*(k)) is that they are dispersion-less (that 
is, completely flat). 

Besides the Andreev bound states within the gap, the 




^ — I 1 1 -u.;j-^ ~ — I ■ 1 

0,00 0.25 0.50 O.OQ 0.25 0.50 



FIG. 5: Spectral function for odd-parity intra-orbital s-wave 
pairing, for model (I) ((a) and (c)) and model (II) ((b) and 
(d)). For (a) and (b) ((c) and (d)), the topological surface 
states are well separated from (merged into) the bulk conduc- 
tion band at the chemical potential. Spectrum along other 
directions (crossing the F point) are qualitatively identical. 
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FIG. 6: Spectral function for odd-parity inter-orbital triplet 
pairing in the equal spin pairing channel, for model (I) and 
A'^' (k)=iAoSo <8> ""2. (a) and (c) are along the kx direction 
while (b) and (d) are along the ky direction. The parameters 
are as shown on the figure. 
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redistribution of spectral weights arising from the particle 
hole mixing is also very interesting. An important fea- 
ture is the appearance of a linear band beyond the Fermi 
momentum and below the chemical potential, existing as 
a particle hole symmetric band of the original topologi- 
cal surface states in the normal phase. Once a bulk su- 
perconducting pairing forms in the topological insulator 
itself (and not in the intercalated copper) in Cua;Bi2Se3, 
such a linear dispersive band is always there, no mat- 
ter a gap opens or not in the topological surface states. 
To see the above feature more clearly, we show in Figs. 
7(a) and 7(b) the energy distribution curves (EDC) for 
several typical wave vectors for two typical pairings and 
parameter sets, which are same as those of Figs. 1(a) and 
3(a), respectively. The linear band mentioned above ap- 
pears as a well defined peak slightly below the chemical 
potential. As the wave vector increases and shifts away 
from the Fermi momentum, the peak deviates linearly 
from the chemical potential and the height and width of 
it both decrease rapidly, which is in agreement with the 
fact that superconducting pairing forms only close to the 
chemical potential. The integrated weight of the linearly 
dispersive peaks in Figs. 7(a) and 7(b) are shown in Fig. 
7(c), as a function of the wave vector. If the gap in the 
superconductors realized from a topological insulator is 
larger than what is reported in Ref. , the above linear 
dispersive structure in the EDC could be detectable by 
ARPES for the wave vectors close enough to the Fermi 
momentum, fiol, ^2d\ Then this well defined peak structure 
arising from the topological surface states could be used 
as a good indicator of the formation of superconducting 
correlation in Bi2Se3 and the involvement of the topolog- 
ical surface states in the superconducting phase. 

From the above results for five different pairing sym- 
metries, we observe a simple rule for the existence of 
nontrivial surface Andreev bound states. For odd-parity 
pairings, when a full gap opens in the continuum part of 
the surface spectrum but no gap opens in the topological 
surface states, a band of surface Andreev bound states 
would arise which traverses the bulk pairing gap. This 
criterion is verified also by calculations for other super- 
conducting pairings not presented here. The results in 
this subsection are summarized in Table I. 



D. the surface Andreev bound states 

For some superconducting pairings, such as the p±ip 
wave pairing, it was known that Majorana fermions ex- 
ist as gapless surface or edge modes. [4Q|] According to an 
argument by Linder et al, all zero energy Andreev bound 
states emerging from the nondegenerate (that is, on each 
surface) topological surface states should be Majorana 
fermions. |4l| As regards our case, in the parameter re- 
gion where the topological surface states are well sepa- 
rated from the bulk conduction band, one observation is 
that the surface Andreev bound states in the gap region 
connect continuously to the topological surface states in- 
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FIG. 7: EDC for three typical wave vectors for (a) even-parity 
intra-orbital singlet pairing (with parameters same as in Fig. 
1(a)) and (b)odd-parity inter-orbital triplet pairing (with pa- 
rameters same as in Fig. 3(a)) within model (I). The param- 
eters are as shown on the figures. (c)The integrated weight 
of the linearly dispersive peaks slightly below the chemical 
potential with parameters corresponding to (a) and (b), re- 
spectively. 



herited from the normal phase (e.g.. Fig. 3(a)). Since the 
topological surface states are spin polarized helical and 
nondegenerate, the surface Andreev bound states on each 
surface should also be nondegenerate. Then according to 
the arguments by Linder et a^(4l|. the zero energy An- 
dreev bound states presented in the above section should 
also be Majorana fermions. 

To see more clearly the properties of the surface An- 
dreev bound states, we perform numerical calculations 
on a finite layer superconducting film. As an example, 
we will analyze the odd-parity inter-orbital triplet pair- 
ing channel described by A(k)=Ao.si ® (T2, within model 
(I). Fig. 8(a) shows the dispersion for a fifty layer film. It 
reproduces all the basic features in the spectral function 
(see Fig. 3(a)). Enlargement of the low energy dispersion 
(Fig. 8(b)) shows that dispersion of the surface Andreev 
bound states is linear close to the F point. For all the 
bulk pairings studied above, the dispersion for the su- 



TABLE I: Summary of results for the bulk pairings considered 
explicitly in the present work. Results for two models, (I) and 
(II), are compared. 'TSS' and 'ABS' are the abbreviations 
for 'topological surface states' and 'Andreev bound states', 
respectively. ' + ' and ' — ' means the even-parity and odd- 
parity pairings. For 'Gap in TSS', 'Y' and 'N' represents that 
a gap could and could not open in the topological surface 
states. For 'ABS', 'Y' and 'N' denotes that Andreev bound 
states exist and do not exist on the surface for a certain bulk 
pairing. 



A(k) 




iS2 ® o-Q 


iS2 ® (Ti 


Si (g) (72 


is2 ® 0-3 


iSo (72 


P 




+ 


+ 








Gap in TSS: 


(I) 


Y 


N 


N 


Y 


N 


Gap in TSS: 


(11) 


Y 


N 


Y 


N 


N 


ABS: (I) 




N 


N 


Y 


N 


Y 


ABS: (II) 




N 


N 


N 


Y 


Y 



perconducting film reproduces well the features of the 
corresponding spectral function. Figure 9(a) shows the 
wave function amplitudes for the surface state localized 
on the top several layers. The corresponding behavior 
for the topological surface states in the normal phase 
is presented in Fig. 9(b). The decay behavior of the 
surface bound states into the bulk in Fig. 9(a) is seen 
to change continuously from oscillatory exponential de- 
cay in the gap region |23| to monotonic exponential decay 
outside the gap region. That is, the state changes from 
particle-hole mixed superconducting quasiparticle to the 
topological surface states in the normal phase. 

In this paper, we have always been discussing a ho- 
mogeneous phase both in the bulk and on the surface. 
However, in the presence of exotic surface excitations as 
vortices, novel Majorana fermion modes may appear in 
the vortex core even for bulk pairings with no surface 
Andreev bound states. There have already been many 
papers focusing on this possibility, which usually start 
from the effective model. 43, 44, 45] 

The Andreev bound states that appear in many differ- 
ent superconducting pairings as shown above confirm the 
idea that superconducting states realized in a topological 
insulator are very probable to have nontrivial topological 
characters. The Andreev bound states, if exist, should 
be easily detectable in a tunneling type experiments as 
a well defined zero energy peak. Another way to detect 
the Majorana fermions as zero energy Andreev bound 
states is to take advantage of the various phase sensitive 
transport devices proposed to produce and manipulate 
the Majorana fermions. [itI l46j 

IV. SUMMARY 

In this paper, we have discussed the surface spectral 
function of superconductors realized from a topological 
insulator, such as the copper-intercalated Bi2Se3. These 
functions are calculated by projecting bulk states to the 
surface for two different models used previously for the 
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FIG. 8: (a) Dispersion of a fifty layer superconductor emerg- 
ing from model (I) and for the odd-parity inter-orbital 
triplet pairing A(k)=AoSi (g) cr2. Parameters used are ^=0.9, 
Ao=0.1, and m=-0.7. (b) An enlargement of the small wave 
vector and low energy part of (a). 



topological insulator. Dependence of the surface spec- 
tra on the symmetry of the bulk pairing order parameter 
are discussed with particular emphasis on the odd-parity 
pairing. When an odd-parity pairing opens a full gap in 
the bulk, but not for the topological surface states, zero 
energy Andreev bound states are shown to appear on the 
surface. When the topological surface states are well sep- 
arated from the bulk conduction band, the redistribution 
of spectral weight induced by the onset of superconduc- 
tivity produces a linearly dispersive peak structure be- 
yond the Fermi momentum and below the chemical po- 
tential. It is proposed as a criterion for confirming that 
superconductivity occurs in the Bi2Se3 (and not in cop- 
per) and the topological surface states are involved in the 
superconducting phase. The zero energy surface Andreev 
bound states are argued to be Majorana fermions. 
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FIG. 9: (a) Decay of the wave function amplitude with layer 
number for the surface Andreev bound states localized on 
the top several layers, results are obtained for the supercon- 
ducting state emerging from model (I) and for the odd-parity 
inter-orbital triplet pairing A(k)=AoSi ® cr2- The parameters 
are /i=0.9, Ao=0.1, and m—-0.7. The Fermi wave vector in 
the kx direction is about O.IOtt. (b) Decay of the wave func- 
tion amplitude with layer number for the topological surface 
states localized on the top several layers, for the normal states 
of model (I) with m—-0.7. The two insets show enlargements 
of the small amplitude regions of the two figures. 
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